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Abstract: We consider dijet production in the region where symmetric cuts on 
the transverse energy, Et, are apphed to the jets. In this region next-to-leading 
order calculations are unreliable and an all-order resummation of soft gluon effects 
is needed, which we carry out. Although, for illustrative purposes, we choose dijets 
produced in deep inelastic scattering, our general ideas apply additionally to dijets 
produced in photoproduction or 77 processes and should be relevant also to the study 
of prompt di-photon Et spectra in association with a recoiling jet, in hadron-hadron 
processes. 
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1. Introduction 

Measurements involving jet production and comparisons of the corresponding rates 
and distributions with QCD calculations have provided some of the best means for 
testing perturbative QCD. As an example, final states involving two or more jets 
have been extensively studied at HERA by the HI and ZEUS collaborations [1,2] 
over a wide kinematic range. At high photon virtualities Q^, comparisons of dijet 
cross-sections and distributions with next-to-leading order (NLO) calculations [3-8] 
have yielded precise measurements of the strong coupling ag [9]. On the other hand 
low dijet measurements have been used as means of obtaining information on 
the gluon distribution xg{x) of the proton, complementary to that obtained from 
structure function scaling violations [2]. However such studies are by no means 
unique to the HERA experiment and dijet production has been actively studied for 
77 coUisions at LEP [10] and pp collisions at the Tevatron [11]. 

A feature common to most experimental analyses on dijets is the presence of 
selection cuts which define the phase space for jet production and are generally 
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meant to ensure that the kinematic region chosen is least affected by theoretical 
uncertainties. In experimental analysis of inclusive dijet final states, one usually 
imposes cuts on the transverse energy Ef of each of the two highest dijets. It was 
observed some time ago by Frixione and Ridolfi [3] that NLO calculations for dijet 
rates break down if symmetric Et cuts Ea, Et2 > -Emin are used on the two highest 
Et jets. The breakdown of the NLO calculation was observed to occur because of 
sensitivity to soft gluon emission, which subsequently led to the region of symmetric 
cuts being dubbed infrared unsafe or unphysical. 

Frixione and Ridolfi suggested the asymmetric cuts Et2 > -Emim > -Emin + 
with A not too small compared to the hard scale of the process, which reduced the 
sensitivity to soft gluon effects and resulted in more reliable NLO predictions. Ex- 
perimental data, on the other hand, can accurately be obtained even in the presence 
of symmetric cuts. For example the cross-section for dijet production in deeply in- 
elastic scattering (DIS) has been experimentally measured over a range of values of 
A, down to the symmetric region A = (see e.g. Ref. [2]). 

A plot of the data for cr(A), the dijet cross-section with a given value of A 
(keeping i?min fixed), versus A can be found in [2] and illustrates the problem clearly. 
It shows that the total rate (t(A) is a monotonically decreasing function of A with 
its maximum value a(0) corresponding to the rate with symmetric cuts, A = 0. This 
is clearly expected on the basis of simple phase space considerations: increasing A 
means that less of the total phase space is available and the dijet rate falls. 

The NLO calculation (the program DISENT [5] was used for this purpose in 
Ref. [2]) on the other hand, agrees with the data for larger A values, but as one 
lowers A there is a turnover of the NLO calculation and the corresponding curve 
starts to fall, whilst the data rises continuously. At A = in particular there is 
a significant difference between the data and the NLO estimate crNLo(O). Hence 
at the very point where the measured cross section is largest there is a maximal 
discrepancy of the NLO result with the data, and in the vicinity of this point, a 
qualitative behaviour different from that indicated by the data. Therefore quite 
clearly, a better understanding is sought of the theoretical hmitations that lead to a 
breakdown of the NLO estimate at small A. 

Moreover the problem discussed above is quite general. It also appears when 
one considers, for example, the hadroproduction of a prompt photon in association 
with a jet. The corresponding fragmentation contribution (when the jet emits a hard 
collinear photon) is important as a background for Higgs searches. Here once again, 
placing symmetric cuts on the final state photon and jet Et values will lead to infrared 
sensitivity of the NLO calculation. Alternatively one can consider prompt diphoton 
production in association with a jet and study the photon pair Et distribution [12]. 
Putting a cut on the recoiling jet Et one can investigate the Et distribution of the 
di-gamma pair. Doing this in NLO QCD, an unphysical discontinuity arises at the 
position of the cut, due to fact that in that region soft gluon emission becomes 
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important. For this paper however we shall continue to use dijet production in 
DIS as our illustrative example and will consider extensions to other processes in 
forthcoming work. 

In the present paper we point out that the unphysical behaviour in the dijet rate 
near A = is due to the presence of large logarithms oi Q/A (where — is the 
photon virtuality) in the slope cr'(A) = da/dA. The logarithms in question arise 
from a veto on real gluon emission above scale A, effective in a certain part of the 
dijet phase space, which causes uncanceled virtual corrections to build up above this 
scale. 

While one will obtain double logarithms (two powers of In Q/ A for every power 
of as) from emissions soft and collinear to the incoming parton, one will obtain single 
logarithms a" In" Q/A from soft gluon emission at large angles.^ These logarithms 
cause the slope calculated at NLO, (T^lo(^)) to change sign becoming positive, at 
small A, and divergent at A = 0. This property of the slope is reflected as a leading 
Aas In^ A term in the NLO computation for the total rate (7(A) at small A. Thus 
while cr(A) has a finite value at A = 0, this value is not correctly given by any fixed 
order of perturbation theory. One needs to first resum the large logarithms in the 
slope cr'iA), to all orders, to obtain a physically meaningful result for cr(A), at small 
A. 

In this paper wc resum soft gluon effects (including hard collinear emission from 
the incoming leg) to all orders in perturbation theory to account for the above large 
logarithms in the slope cr'(A). Our resummation will be in the space of a Fourier 
variable b conjugate to A and we shall resum logarithms in b that, at large b, reflect 
the singular behaviour at small A. These logarithms shall be resummed into a form 
factor E, which can be expressed as 



with L = \n{bQ) and gi and g2 being the leading and next-to-leading logarithmic 
functions. Wc shall refer to this as next-to-leading logarithmic (NLL) or single- 
logarithmic (SL) accuracy. 

Another factor we have to consider however is the conservation of transverse 
momentum. Vectorial cancellations between harder emissions are another way of 
obtaining a small Et difference between the final state jets and this effect also impacts 
the slope a' (A) at very small A. The full answer will be given by the convolution of 
the form factor E with an oscillatory function: 



^In several common jet algorithms, e.g. cone algorithms and their variants, and the inclusive ki 
algorithm [13, 14], there will be no soft and collinear double logarithmic (DL) contributions from 
the outgoing hard partons, provided one recombines partons into jets appropriately, which we shall 
discuss shortly. 



lnS(6) = Lgi{asL) + g2{asL) 



(1.1) 




(1.2) 
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where the sine function accounts for vectorial transverse momentum conservation. 
Once this convolution is performed to obtain the resummed slope one finds that the 
unphysical behaviour is replaced by a smooth behaviour in the limit A — > 0. In fact 
instead of diverging the slope remains finite (and negative as is physically required) 
and for A/Q <^ 1 is proportional to A, i.e. is linear. 

To obtain the maximal possible accuracy, one has to match the resummation per- 
formed here, with fixed order computations that account for subleading logarithms 
and finite corrections (constant pieces and non-logarithmic terms in A). These would 
start at NLO and will be important to get a good description at larger A values. We 
shall postpone dealing with the issue of matching to subsequent work. 

In all of the above considerations, the definition of jets will naturally have a 
significant impact on the answer. One has to choose both a jet algorithm and a 
recombination scheme which details how the kinematic properties of the jet (such as 
its Et) relate to those of its partonic constituents. The results we present here are 
based on the use of a cone-type algorithm also used previously in theoretical studies 
involving dijets [15]. We shall mention some details of this subsequently. We require 
also that particles are clustered into jets using a four- vector recombination scheme 
Pjet = Siejct Pi where p labels four-momenta and the sum runs over all partons that 
constitute the jet. With this recombination scheme the leading logarithmic function 
Qi is independent of the details of the jet algorithm, since it arises purely from 
emissions coUinear to the incoming parton, which will not be recombined with the 
outgoing jets ^ Hence the specifics of the jet algorithm enter the function g2, i.e. at 
NLL level. 

The outline of this paper is as follows. In the next section we shall consider the 
situation at leading (Born) order and introduce all the quantities involved. In the 
following section we illustrate the origin of the soft gluon problem at NLO in more 
detail and calculate the DL divergence at NLO. Next wc shall perform the all-order 
resummation to NLL accuracy, in b space. Subsequent to this we shall present our 
final numerical results and a discussion illustrating the main features thereof. Lastly 
we shall conclude while mentioning some planned future developments and work in 
progress. 

2. Dijet production at leading order 

Let us consider the production of two final state jets in the Brcit frame of DIS. To 
leading order, the dijets are just two partons labeled k and r (see figure 1). We also 
denote with p the incoming parton four-momentum and with q the four-momentum 

^What we actually need, for our calculations to directly apply, is a recombination scheme that 
vectorially adds the three-momenta of partons within a jet. Then the jet Et is just defined as the 
magnitude of the corresponding jet transverse- momentum vector Ptjet- 
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of the virtual photon. Further one imposes the asymmetric cuts 

En > -Emin + A , Et2 > -Emin , (2.1) 

where Eti denotes the transverse energy of the i*'* jet (and to this order parton r or 
k) with respect to the photon axis in the Breit frame. 




Figure 1: The Born kinematics of dijet production in DIS. 

The general expression for the total rate (with the above cuts imposed) for dijet 
production in x space can be written at leading order as 

cToix, Q, E^^, A) = £ ^Co(e, Q, E^^, A)/ (^|, Q'^ , (2.2) 

where the subscript '0' denotes the fact that the above result is at Born level, x 
denotes the usual DIS Bjorken variable while ^ = Q^/ (2(p ■ q)) is the fraction of the 
incoming parton momentum p carried by the struck parton. In the above formula 
we have used Cq to denote a generic coefficient function, implicitly including in it 
the transverse and longitudinal components (for simplicity we confine the discussion 
to virtual photon exchange only). Accordingly f{x/^) has been used to denote the 
parton density and includes the dependence on parton flavour and charge. We shall 
always consider the renormalisation scale and the factorisation scale as equal to Q, 
but variations around this value can be trivially accounted for. 
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The leading order coefficient function Co can be obtained by integrating the 
squared matrix element over the desired phase space as below: 

X Q{En - {E^n + A))e{Et2 - E^i^) . (2.3) 

Here M^{p,q,r,k) is the leading order matrix element squared for the hard dijct 
production at leading order and is made up of invariants constructed from the various 
parton momenta. It differs according to whether the subprocess we consider involves 
an incoming quark or gluon but the general form above applies in both cases. We 
shall avoid displaying this dependence explicitly as well as the dependence on i?min 
and Q, in what follows below. In (2.3) Tq and ko denote the final-state particle 
energies. 

Integrating away various components, using the delta function and noting that 
at Born level \ft\ = rt = \kt\ = En = Et2, we are left with 

where and k^ are the corresponding parton longitudinal momentum components 
along the photon axis, now fixed in terms of the components of ft, and we identified 
each outgoing parton with a jet. 

We now introduce the slope (t'(A) = da/dA. At leading order this is just 

a'oix, Q, E^^, A) = £ ^C;(e, Q, E^^, A)/ (^|, Q'^ , (2.5) 

with the coefficient function obtained by differentiation of (2.3) with respect to A, 

8C r fl'^f 1 

A) ^ ^ ^ - y ^^^-^—^^Hr. - (E„,„ + A)) (2.6) 

This integral can be performed (with any additional cuts such as one on the 
interjet rapidity) and has a finite value as A — > 0. Moreover at this order the slope 
is negative at all A as one requires on physical grounds. As we shall see this is no 
longer the case at NLO. 

3. Soft gluon effects at NLO 

The aim of this section will be to discuss the kinematical constraint on soft gluon 
emission, that arises in the region of small A/Q, when one moves beyond the leading 
order eq. (2.5). This constraint results in logarithmic enhancements and we shall 
exphcitly compute the DL behaviour ctsln^Q/A, that first arises at NLO. Before 
that we discuss the relevant kinematics. 
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3.1 Kinematics 



Moving to NLO we have to treat additional gluon emission. If the additional gluon 

is rccombincd with an outgoing hard parton by a jet algorithm and with four vector 
addition, it docs not cause a mismatch in the jet transverse energies, En = Et2, and 
the jets are exactly back-to-back in azimuth. In this region the soft gluon contri- 
bution cancels with virtual corrections. If however the gluon is not recombined into 
an outgoing jet, e.g. it is near the beam (incoming parton) direction, it contributes 
to a transverse energy mismatch between the jets. Now the soft gluon effects do not 
cancel fully with virtual corrections and large logarithms appear. 

To explicate this, we write the four-momenta of the outgoing partons as (here 
we explicitly identify partons r and k with jets and assume k' is not recombined with 
them) : 

n = Etiil,0), (3.1) 
h = Etaicos^n ± e) , sm{TT ± e)) , (3.2) 
k't = k'f{cos4>, sin (f)) , (3.3) 

where the two jets are almost back-to-back in the transverse plane, since if k' is soft 
the hard parton recoil e is small. Prom transverse momentum conservation one gets 

= £2 + E^^ - 2EaEt2 + e^E^^E.^ + Oie^EnE^) ■ (3.4) 

Additionally using = [k'^ / Ef2) sm^ (p (assuming En Et2 in the soft emission 
limit, discarding subleading terms and terms that are important over only a para- 
metrically small interval in (f) one has simply that 

k',\cos^\^%\^\Ea-Et2\. (3.5) 

The terms we neglected will not contribute at the NLL accuracy we aim for 
in this article. Thus at our level of accuracy, the mismatch in jet Et arises from a 
particular component of soft gluon transverse momenta (the component along the jet 
axis in the plane transverse to the photon axis). Now that we have established how 
precisely soft gluons flowing outside the jets contribute to an Et mismatch between 
them, we can consider what happens due to the placement of Et cuts on the high E^ 
dijets. 

In this regard it is helpful to examine the diagram in figure 2, which depicts the 
situation in the E^., Et2 plane. The shaded rectangular region shown, corresponds to 
Eti > ii^min+ A, Et2 > -Emin and hcucc dcuotcs the region allowed by the experimental 
cuts. Let us consider the contribution from points along the dotted hue En — 
Emin + A, which contributes to the slope as is evident from (2.6). At Born level 
additionally, we are kinematically constrained to be on the line En — Ef2 (the lower 
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Figure 2: Phase space diagram in the Eti,Et2 plane. 

of the diagonal lines in the figure) and hence the contribution we obtain is from the 
single point En = Et2 = E^in + A. 

Now consider the emission of a soft gluon k' that causes a small mismatch En = 
Et2 + S, depicted by the upper of the diagonal straight lines in the plot. For this soft 
gluon to contribute, it must intersect the dashed horizontal line En = E^^^+A, inside 
the allowed region (and hence must pass through the shaded triangular region). If the 
soft gluon has energy (more precisely as we discussed the component kx) 6 > A then 
it is pushed outside the allowed region and hence vetoed. A veto on soft emissions 
above some small scale causes uncanceled virtual contributions at that scale, which 
in turn leads to logarithmic behaviour in A/Q. This behaviour will be manifest in 
the derivative of the total rate with respect to A, which receives contributions only 
from points on the dotted line En = -Emin + A. 

Such a DL behaviour in the slope of the dijet rate (plotted against A) is present 
in the fixed order computations with DISENT [5] and causes an unphysical turnover 
of the theoretical calculation as we move to the small A region. A resummation of 
soft gluon effects is thus required to restore the correct physical behaviour seen for 
example in the ZEUS data [2]. 

3.2 Soft one- loop calculation 

We shall now compute explicitly the above described soft gluon behaviour at NLO 
(one-loop). Adding a real soft gluon with four momentum k' (with energy uj' ^ Q) to 
the Born system (which in both incoming quark and gluon channels is a configuration 
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with three hard partons, two fermions and a gluon), we have the real contribution 
to the NLO coefficient function^ 



2(2.)3.„ 2(2.)3A, 2(2.)3n ^ ^ - ^ — '•^''^ 

X s[k']e(E„ - (£;„,„ + A))e(B,2 - £;„,„). (3.6) 

Prom now on, to ease presentation, we shall only consider the calculation for the 
incoming quark channel with partons r and k being respectively an outgoing quark 
and gluon, exactly as depicted in figure 1. For the final results we include all possible 
configurations, i. e. the contribution from the incoming gluon channel as well. 

S[k'] is the three particle {qqg) antenna pattern, which (for our chosen channel) 
is given by: 



Wrk + Wpk - -TJ^W. 



(3.7) 



where 

with Nc = Ca , and being the strong coupling such that 5'^/47r = a^. 

At NLO the leading real soft gluon contribution to the slope a' can be obtained 
by differentiating with respect to A in eq. (3.6), which then gives 



8{2ny J Ir-sko-ksrol k'^ 

X 5^{kt + ft + k't)e{Ea - E^,^)5{Ea - {E^m + A)) (3.9) 

Using vectorial kf conservation, we can remove the delta function recalling that 
it leads to the condition \kx\ — \Eti — £'42 1- Prom the region En > Et2 and after 
introducing the virtual emission of k' (which has weight —S[k']), we obtain the 
following expression, vahd at small A (see figure 2): 



X 



^s[k'mA-\k'j)-i] , 



(3.10) 



where the integration in k' extends to the region where k' is not recombined with k 
or r to form a jet. Note the step function that restricts \k'J in the real piece but not 
in the negative virtual contribution. 

^Wc have neglected the dependence of the parton distributions on the soft gluon k' which we 
are allowed to do for soft emissions . We shall correct for this effect when including next-to-leading 
(single) logarithms arising from hard emissions coUinear to the initial state partons, where one 
cannot make this simplification. 
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On combining with the leading (Born order) result and wishing to retain only 
logarithmic terms in A, the slope of the total rate at small A can be described at 
NLO as 

C;Lo(e, A) = + = A)iyNLo(A). (3.11) 

It is simple to calculate Wt<si^o{A), although the full answer will depend on the 
jet-algorithm employed. However the leading DL term in Q/A follows by collecting 
the coUinear singularities along the incoming leg p in S[k'], assuming an algorithm 
that ensures that gluons coUinear to the final state hard partons do not contribute, 
and we obtain (recall that we are considering the incoming quark channel) 

W^^o{A)^l + J ^^0^5[fc'][e(A-|A;;|)-l]«.l-C^^ln2 !+•••, (3.12) 

where the dots denote SL terms that we shall compute later, and constant pieces 
and terms that vanish at small A, which need to be accounted for by matching. 

To summarise, at small A one inhibits the radiation of real soft gluons and 
this gives rise to the DL behaviour in the slope a' via the double logarithm in 
the coefficient function C^lo- This DL behaviour is responsible for an unphysical 
turn over of the NLO calculation at small A [2] . In order to cure this pathological 
behaviour one has to perform a resummation of such soft gluon effects to all orders. 
This will be the aim of the next section. 



4. All order result 

The extension of the NLO result cq. (3.11) to all orders has two main elements: the 
computation of multiple all-order soft giuon emission from a three particle antenna 
comprising two final state and one incoming hard parton (global component) and 
accounting for correlated emission outside the final state jets from soft gluons within 
the jets (non-global logarithms) [16-18]. At NLL level the details of the jet algorithm 
become relevant, so that before proceeding with the calculation it is useful to have 
a procedure to refer to. One can, for example, use the cone algorithm introduced 
in [15], which samples the phase space for sets of particles flowing into cones of fixed 
angular size, 6, and is defined as follows: 

1. given a set of final state momenta X — {qi}, for any subset x oi X compute 
the unit vector 

^ ; (4.1) 

2. consider the set x' of all particles that flow inside a cone of half- angle S centered 
around fix; 

3. a jet is any set of particles x for which x' — x. 
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With such a procedure one can, of course, generate any number of jets. However 
we are not concerned with the detailed jet structure of the final state, for example 
identifying the number of jets. We are just concerned with the two highest jets 
generated by this procedure and the Ef mismatch between them. We shall consider 
S small compared to 1, but not so small that one needs to resum logarithms in the 
cone-size, cts In 1/5 <S 1. Equipped with this procedure we turn to computing the all- 
order resummed result. We begin by examining the multiple independent emission 
(global) contribution. 

4.1 Global component 

To derive this part of the result, one just considers that multiple soft gluons are 
emitted independently of one-another, each following the antenna pattern S[k']. 
Secondary parton splitting is built-in via the running coupling. However, in the 
present case, this approximation misses a subset of the next-to-leading logarithms, 
which arise from correlated as opposed to independent emission (non-global loga- 
rithms) [16, 17]. In other words gluons emitted outside but near the jet boundaries, 
feel the dynamical influence of relatively harder emissions that are inside the jets 
and this configuration generates next-to-leading logarithms. 

The global part of the answer has both leading (double) and next-to-leading (sin- 
gle) logarithms. The leading logarithms shall arise from emission soft and collinear 
to the incoming hard parton. The next-to-leading logarithms in the global piece will 
arise from two sources. Firstly soft radiation at large angles to the three hard emit- 
ters is a source of single logarithms, which will contain a characteristic dependence 
on the geometry of the three-pronged hard antenna. An additional source of single 
logarithms is hard emissions quasi-coUinear to the incoming parton leg, which we 
shall also include. 

To derive the all orders global result we start with the observation that the 
function W{A) can be defined by extending the NLO result eq. (3.12) to all orders, 
using the factorisation property of multiple soft gluon ensembles: 



where the above integral includes contributions only from gluons i ^ 1,2, i.e. that 
fly outside the final state jets 1 and 2. 

We first factorise the step function as follows: 




-S[k[]) e(A-|^fcV|)-l , (4.2) 




1^1,2 




(4.3) 
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where we used the Fourier transform of the delta function 

5{x) = / — exp[-ibx] . (4.4) 
J-oo 27r 

This allows to simplify the above further to read (valid near A = 0) 

C'it A) = Cl, + C[ + ---^ ait A)W{A) , (4.5) 

with 

1 f°° dh 

W(A) = - ^sin(6A)exp[-i?(6)]. (4.6) 
Here the 'radiator' R(b) is given by 

^^^^ ^ J 2(2^'^[^'^ " ^M^bK]) ■ (4.7) 
To NLL accuracy we can simplify the radiator via the replacements 

1 - exp[ibk'^] 1 - cos(6A;;) ^ Q {K\ cos (/)| - 1/6) , b^be"'^ . (4.8) 
This allows us to achieve our final form for VF(A), which reads 

WiA)^- sin(6A)exp[-m)]. (4.9) 

Jo 

with radiator at NLL accuracy given by: 

/d^k' 
^^^^s[k']e{K\cos<p\-i/b) . (4.10) 

We emphasise that the integral over k' is in the region where it is not recombined 
with an outgoing jet. 

We now proceed to the computation of the radiator in the particular case of an 
incoming quark. During the calculation we will discuss how the results obtained can 
be generalised to the incoming gluon case. 

4.1.1 Leading logarithmic result 

At the leading logarithmic (LL) level the radiator is easy to compute since one has 
to consider just radiation collinear to the initial state parton. Radiation collinear 
to either r ot k will be clustered into a jet and hence in a cone-type or inclusive 
kt algorithm the only source of leading logarithms will be from this initial state 
radiation. 

Collecting the collinear singularities along the incoming direction in S[k'] and 
using eq. (4.10) we arrive at 

i?BL(6) = 2C, I In (I) g e(^;| COS 01 - b-') = In' {bQ) + ■■■ 

(4.11) 
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where to extract the double logarithms we froze the coupling at scale Q. 

Doing the above integral with the running coupling converts the DL contribution 
into a LL function given by 

^20j'§A'l. ^ -r/J^ , (4.12) 

Ji/b K K din 6 

where we have performed an expansion of Rhhib) around b, which we are allowed 
to do since what is left is a subleading contribution. At NLL accuracy the radiator 
has the following expression: 

i?LL(&) = i^^i(A) + /2(A), (4.13) 

where gi is the leading logarithmic result, /2 is a piece of the NLL contribution g2 
and L = ln{bQ) while A — Poas{Q)L. The leading logarithmic result gi reads 

^^^^^^"2^t2A + ln(l-2A)] , (4.14) 

while /2 is given by 

,,,, KCf 2A \ /3iCf /I ln(l-2A) + 2A^ 

(4.15) 

where and Pi are the first two coefficients of the QCD beta function: 

_ IIC^ - 2nf _ \1C\ - bCAUf - SC^Uf 

12;^' 24;^^ ■ ^^-^^^ 

In order to account for all NLL contributions coming from soft and collinear radiation, 
the coupling constant in the k[ integral in (4.12) has to be taken in the physical CMW 
scheme [19], which is related to the MS scheme by the relation 

The logarithmic derivative of Rhhib) in (4.12) can be obtained by differentiating only 
the gi piece of Rhhib), since what is left is subleading. 

Note that although we have labeled the piece of the radiator computed here as 
-Rll, we also include in it NLL terms arising from the running coupling and change 
of scheme. It is perhaps better to think of this as a DL piece (arising from soft and 
collinear emission, while the next-to-leading logarithms we compute subsequently 
are either pure soft or pure coUinear SL (q;"L") effects. 
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4.1.2 NLL soft contribution 

Having computed the leading logaritlimic piece of our answer, which is independent 
of the details of the jet definition {e.g. the cone size) we now turn our attention to 
the NLL terms arising from the independent emission (global) part of the answer. 
Non-global NLL terms to do with correlated emission will be treated in the next 
subsection. To include global NLL terms we need to compute eq. (4.10) using the full 
form of the soft function S[k'], rather than just collecting the coUinear singularities 
on the incoming leg, as we did for the leading logarithmic terms. This will enable us 
to correct for SL terms arising from soft, coherent interjet radiation. Additionally we 
have to treat the dependence of the variable on the azimuthal angle 0, which could 
be discarded at leading logarithmic level, and correct for hard collincar emission. 

We first treat the (f) dependence. To NLL accuracy eq. (4.10) can be written 
as (performing a Taylor expansion about | cos0| = 1 of the full result and retaining 
terms only up to NLL accuracy) 

where the ellipsis denotes terms beyond NLL accuracy which would be produced by 
taking higher derivatives or taking the derivative of any piece of R{b) that is not 
leading logarithmic. The function -Rll(&) was computed in eq. (4.13) (in fact the 
only relevant contribution to the derivative above will be from the gi piece of it^LL 
and we shall discard other terms) while J^"^ \n{ \ cos0|)|^ = — In 2. 

This leaves the first term on the right hand side of the above equation which con- 
tains both the already computed leading-log terms and next-to-leading logarithms 
yet to be accounted for. In order to compute it fully one can treat each dipole in 
S[k'] in turn. For example let us consider the rk dipole where r and k initiate the 
final state jets. We shall take r and k to be an outgoing quark and gluon respectively 
(see figure 1) and the corresponding contribution to the first term on the right hand 
side of (4.18) is 

where [rk] (for instance) denotes 1 — cos6rk, and the scale (k^)^ = 2{rk'){kk')/{rk) 
is the transverse momentum (squared) of the soft emission k' with respect to the 
emitting dipole pair. This scale naturally emerges when one considers the collincar 
splitting of gluon k' into two offspring partons with similar energies [21]. Note that 
the angular integration over the directions, nj,, of k' in the above equation is con- 
strained such that k' is outside a cone of angular size 5 around the outgoing hard 
partons r and k, k' ^ 5r,k- This is of course an approximation since, in the chosen 
algorithm, 5 is really the allowed opening angle wrt the energy weighted centroid 
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of the outgoing hard parton and the emitted soft parton. The correction terms are 
proportional to powers of the gluon energy and we are entitled to neglect them here, 
for the small A behaviour. 

In principle the integral in eq. (4.19) is rather cumbersome to evaluate in full de- 
tail. However one can simplify the situation to extract the LL and NLL dependence. 
In particular the rk dipole does not contribute any leading logarithms since these 
arise from emissions that are soft and coUinear to the incoming leg p. Hence the 
contribution from the rk dipole is at most NLL in b (recall that we do not attempt 
to resum logarithms in the cone size S). The single logarithms in question arise from 
the pole in the integration over energy uj' and are a soft wide-angle contribution. To 
extract this piece we can simplify (4.19) to give 

R.M^cJ^^e(u,'-l/m^?^, (4.20) 

where we set iVc = Ca- Notice that we felt free to mistreat the argument of ag 
but retained its essential dependence on the energy uj' in doing which we neglected 
a constant of proportionality which would produce only NNLL terms beyond our 
accuracy. In writing the above we also made use of the result 

and discarded terms involving the ratio of the cone-size (squared) to the interjet 
separation, which we shall do throughout this paper. However note that we can 
retain the full dependence on cone size, in this wide-angle global piece, by employing 
the exact formula mentioned above rather than retaining simply the logarithmic 
dependence on cone-size S. Note also the dependence on the geometry [rk] of the 
underlying dipole emitters, that is typical of soft interjet radiation [21]. 

Similarly one can compute the other dipoles pr and pk. Here we will also en- 
counter leading logarithms from when k' is near the incoming leg p, and in this region 
the argument of the running coupling will reduce to k^ rather cu'. However we have 
already computed these leading/double logarithms in (4.12) and the remaining NLL 
piece of the answer will once again be obtained by arguments along the lines above. 
Assembling the contribution from all dipoles we have below the full soft contribution 
to the radiator which can be expressed as 

/ 2(0^*^^^]® + ^^^^^^^ ' ^^-^^^ 

with i?LL as given in eq. (4.12) and -Rnll being the soft global NLL contribution to 
the radiator: 



R 



■NLL 



2C..(?M).C..(™),,., 



(2t). (4.23) 
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We introduced above the SL evolution variable t: 



A/b u 2n inpo 1 - 2A ^ ^ 

Recall that we treated p and r as quarks and A; as a gluon, but for the final results 
we sum over all configurations with appropriate modifications to the above form. 

Notice that the SL contribution (4.23) depends not only on the geometry of the 
three-parton antenna (specifically on the angles between hard emitting partons) and 
the cone size S, but also on the incoming parton energy Ep. The additional term 
Cp^'O-Ep/Q'^ accounts for the fact that a soft gluon collinear to p has energy up to 
Ep — Q/{2^), and not Q, as one would infer from (4.12). 

There are still two sources of single logarithms missing from the above answer. 
The first source of single logarithms is from non-soft emissions almost collinear to 
the incoming parton p, which we shall include next, to complete the global piece of 
the calculation. The other piece we will need is the non-global term arising from soft 
correlated emission of k' from gluons included within the jets. 

4.1.3 NLL terms from hard collinear emissions 

Here we note that multiple hard emissions on the incoming leg p also contribute a 
class of single logarithms, precisely as in the case of DIS event shape variables [22-24]. 

To exponentiate this piece one has to turn to Mellin space with respect to the Bjorken 
X variable. However we shall directly note (referring the interested reader to the 
manipulations described for instance in [22]) that restricting the kt of hard collinear 
emissions on the incoming leg, one essentially restricts the DGLAP evolution of the 
structure function to the scale 1/6^, rather than Q^, f{x/^, Q^) — > f{x/^, 1/6^). One 
also needs to change the scale in the calculation for i?LL such that the corresponding 
virtual corrections are properly treated. This leads to the replacement of Q in the 
integrand of eq. (4.12) by the factor Qe"^/^. Thus an additional term appears in the 
radiator which is given by 

and the remaining piece of the hard collinear emission is embodied in a change of 
scale of the parton densities as mentioned before. In the incoming gluon case the 
hard collinear contribution can be obtained from (4.25) by simply replacing —3/4 
with —ttPo/Ca, with Pq defined in (4.16). 

4.2 non-global component 

So far we accounted only for independent multiple soft gluon emission, with correc- 
tions for hard-coUinear emission. For a class of observables that typically involve 
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angular cuts in the phase space, the independent emission approximation is not suf- 
ficient to generate the full single logarithms [16, 17]. Our observable is indeed such 
a non-global observable, due to the fact that it is sensitive to radiation only outside 
the jets. Hence a soft emission near the jet boundary, which contributes to our ob- 
servable, can resolve relatively harder emission inside the jets, i.e. the jet structure, 
at NLL level. 

We now account for the non- global piece of the final result, for the slope in 
A of the dijet rate (7(A). As we said, non-global single logarithms, Q;"ln"(6(5), 
arise from soft emissions that fiy inside the jets defined by cones, which themselves 
emit outside the jet region. On a heuristic level, emission from jets at large angles 
compared to the angular extent of the jets themselves, will see only the total colour 
charge of the system of partonic emitters that constitute the jet. Thus it follows from 
coherence properties of QCD radiation that when one considers the relatively small 
cone approximation (jet cones significantly smaller than the interjet separation) the 
non-global component will arise separately from each cone boundary. We illustrate 
our ideas by first performing below an analytical computation of the leading In^ bQ 
non global piece and follow that by considering non-global effects at all orders. 

Before we move on we should however mention that non global logs are rather 
sensitive to the exact details of the jet algorithm employed. For instance in our case 
(cone algorithm), a situation could arise where a soft parton may form a jet with 
the hardest (high Et) parton and also with a softer parton itself outside the high 
jet. The decision on how to attribute the common energy between the jets will affect 
the size of the non global component. In what follows we ignore this complication 
and stick to our previous definition of the high E^ dijet, which will mean we take a 
scenario where non-global logs make a maximal impact. 

In order to proceed we need to consider the emission of two soft gluons by a 
hard three-particle antenna. This has a colour dipole structure identical to that of 
the single emission term eq. (3.7) except that for each dipole term Wab of eq. (3.7) 

(2) 

one inserts the result for emission of a soft two-parton system by the dipole ab, w\^l , 
that is also the relevant function in the two-jet case (see [21,25,26]). In fact for 
the non-global term we shall need to consider only a specific piece of the correlated 
two-parton emission term, corresponding to the situation when the two soft emitted 
partons are energy ordered, uji 3> a;2. Its detailed form will be mentioned below. 

We thus again consider a generic dipole ab formed by two of the three hard 
partons that are present in the underlying hard event. We parametrise the four- 
momenta of these hard partons as below, along with those of the emitted two soft 
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gluons ki and k2: 



p„ = 0,0,1), 

Pb = Eb{l, 0, sin 9ab, cos Oab) , ^g-j 
fci = ci;i(l, sin 9i sin 0i, sin^i cos 0i, cos^i) , 
^2 = a;2(l, sin ^2 sin 02, sin ^2 cos 02, cos ^2) , 

and assume strong energy ordering uii ^ UJ2 as is required to generate the non-global 
piece. To trigger the non-global contribution, the harder parton ki flies inside the 
jet cones while k2 is outside. For the dipole ab we assume that parton b is the hard 
incoming parton p while a is an outgoing hard parton (r or k), which initiates a jet 
with angular size S. In general both legs of the dipole ab can correspond to the two 
outgoing jets {e.g. a = r,b = k) and this contribution will be treated subsequently. 
For the present case however non-global logs are generated by the configuration 

cos 6*1 > cos 5 , cos ^2 < cos (5. (4.27) 

The leading non-global contribution is then obtained similarly as in [16, 17] by 
considering the integral (once again we only need to consider the dependence on the 
energy cu, since this is a pure soft piece) 

V27ry Jo Jo ^2 27r 27r H ^ ' ' ' ' ^ ' 

and we denoted by w^^^,] angular dependence of w^^^ . Performing the energy 
integrals is trivial and gives — (^)^ ln^(6(5). To work out the coefficient of this SL 
term, we need to integrate w^^^] allowed directions of ki and A;2, rfi and ri2. 

For this we need just the angular dependence of the piece of the correlated emission 
term w^^-* which is given by [25] 

, (2) ( [ab] [ab] [abf \ 

["^] '^VM[12][26] [a2][21][26] [al][16][a2][26]; ' ^ ' 



with, as before, [ij] = 1 — cos 6ij 

lb] 



We first perform an azimuthal average of wPy and get^ 



/ (2) \ 2Ca (1 - cosggfc) e(cosg2 - COSggj,) 

\^[ab]) (1 _ cos 62) (cos 01 - cos 62) (cOS 0i - COS 0ab) ' 

In doing the calculation we have exploited the fact that the only contribution to 
non-global logs arises when the harder gluon ki is emitted inside the jet-cone around 



^We are free to do this since one can discard at the SL level the observable's dependence on 
azimuth and just consider its energy dependence as indicated above. 
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a and the softer gluon ^2 is emitted outside. Note that both gluons are soft compared 
to the hard scale Q. 

Now we need to integrate over polar angles, which gives 




cos Oab 



1 — COS 62 



d cos 62 




1 



d cos 9 



(COS^I — cos ^2) (cos ^1 — COS^ab) 



2Ca(1-cosM 



- Ca-, (4.31) 



independent of the cone size 6. The fact that the result does not depend on the cone 
size illustrates the fact that the contribution from emission of k2 is an edge effect 
arising from the vicinity of the cone boundary and the geometry of the interdipole 
region becomes unimportant, since it corresponds essentially to an infinite interval 
in rapidity between the cone boundary and the other (incoming) emitter (see [17]). 

This is no longer true when one considers the dipole formed by the two outgoing 
partons and in this case there are contributions from each cone boundary. The final 
result depends both on the cone-size 5 as well as the interdipole separation [ah]. 
However if one considers the relatively small cone limit such that one neglects terms 
that vary as 5"^ /[ah] -C 1, the result from such a dipole is simply twice that in (4.31). 
Recall that such finite cone-size effects also arose as corrections to the In 1/5 piece 
of the global single logarithms and we choose to neglect them. 

The final result for the leading non-global term produced by the three-jet sys- 
tem then is given by adding up the contributions from each of the three hard emit- 
ting dipoles with the appropriate colour factors {Nc/2 for a quark- gluon dipole and 
— l/2Nc for a quark-antiquark dipole). Denoting the entire non-global contribution 
by the series 



n=2 

with t the SL evolution variable introduced earlier in eq. (4.24), we have computed 
the first term 5*2 which reads 



Here Ci and C2 are the charges of the partons that initiate the outgoing jets in a 
given underlying hard configuration. For our chosen channel with incoming quark, 
we have Ci + C2 — Cp + Ca- The fixed order result above is correct up to terms 
5"^ /[rk] where r and k are the outgoing hard partons. 

To generalise the above result to all orders one has to consider configurations of 
several wide-angle soft gluons inside the jet cones, that coherently emit a single softest 
gluon outside the jets. At present this can only be done in the large N(. limit which 
reduces the problem to planar graphs, and the corrections to these are suppressed 
by 1/A^c- Therefore one expects these neglected terms to make a difference at the 
10% level to our non-global result [16-18]. 




(4.32) 



-^2 = - 



Ca{C^ + C2)—. 



(4.33) 
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In this limit one can consider that the contribution of each initiating dipole will 
be modified by a non-global factor^ 

e-«ai. ^ e-^'^'Sah (4.34) 

where e~^'^^ is just the independent emission radiator computed earlier, for dipole 
ah and Sab is the accompanying non-global factor (see arguments in e.g. [18]). For 
dijets with opening angles 5 that are small compared to the interdipole opening 
angle [afe] , the contribution Sab coming from each cone boundary will in fact be the 
same as that obtained for the case of a emission into a large (effectively infinite) 
rapidity slice, already derived for the two-jet (single initiating dipole) case [17]. This 
property shows up in the fixed order computation but it was also formally shown 
at all-orders in [18] that for problems with limited energy flow everywhere except 
in small cones around the leading hard partons, the non-global contribution arises 
separately from each cone boundary. This ensures that one can extend the results 
derived previously [16, 17] for emission into an infinite rapidity region, to the present 
case. The corrections to this result will vary as the square of the cone-size 5. 

Putting together the contribution of all the three initiating dipoles (the all order 
non-global result is computed in the large Nc limit) corrects the SL global result 
described earlier: 

g-RW ^ ^~Rib)s{h) , (4.35) 
where S{h) can be parametrised as below [16]: 

S.exp(-(a + C,)C.^(i±(g)*^). (4.36) 

with t as in (4.24) and 

a = 0.85Ca , h = Q.86Ca , c = 1.33 . (4.37) 

We have retained through the above parametrisation the correct Ca{Ci + C2) 
colour structure for the leading term 5*2 calculated earlier, but beyond this leading 
term the result is correct only in the large limit. The paramtrisation above is 
valid for t < 0.7, which is more than sufficient for our purposes. This is because, 
in practice, the b integral receives vanishingly small contributions near and beyond 
this point, which corresponds to PqC^sL = 0.497, very close to the Landau pole value 
of 0.5. In practice however, non-global logs make a very small contribution to the 
overall result. This is due to the fact that they start at 0{al), relative to the Born 
term, and in the region where they may be expected to formally be important (at 
very small A), the b space integral is dominated by the DL behaviour. 

^Strictly speaking this will be true of only the quark-gluon dipoles which survive the large A*",; 
limit. However doing the same also for the large A^c suppressed quark-quark dipole will ensure 
compatibility with the leading order result. Any differences with the correct full answer will of 
course be suppressed as l/N^- 
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5. Results and general properties 

Here we shall present our final result and illustrate its main properties. We can 
express our resummed result in the form (valid at small A) 

a'(A)«£|ae,A)iy(^A,0 , (5.1) 

where we suppressed the dependence on Q and £'min and Cq is the Born level coeffi- 
cient function for the slope defined in eq. (2.6). We redefined the all orders quantity 
H^(A), to include the dependence on the parton densities, 

W (a, ^)=lfj sin(6A)e-^(^)5(6)/ (| 1) . (5.2) 

Here R(b) = i?LL(^) + -Rnll(^) + -^nll(^)' ^^'^ radiator computed in the text, 
in terms of the separate pieces indicated. S{b) is the non-global contribution. 

The most important feature of our answer is the absence of a Sudakov peak when 
one goes from b space to A space. This feature emerges at the DL level itself, and 
hence while SL effects make a significant numerical difference to the final answers, 
the properties pointed out in the following discussion are essentially unaffected by 
neglecting the single logarithms. Therefore for this discussion, the relevant quantity 
we examine is 

W^dl(A) = - / ^ sin(6A)e-^°^('') , (5.3) 
^ Jo b 

where DL indicates that we kept only the DL terms in b space and use for i?DL the 
expression given in eq. (4.11). 

As was pointed out in Ref. [20] the integral above has two distinct regimes. In 
order to study these separately we divide the integral above as follows 

nOO nbo nOO 

/ db= db+ db, (5.4) 

Jo Jo Jbo 

where bo is taken as 1/ A and A = Ae'''^ . The discussion that follows holds also for 
other possible choices of bo, as long as its value is of order 1/A. 

First we shall deal with the second term on the right-hand side of the above 
equation, where the integral is dominated by its lower bound bo- For this term one 
can invert the b transform retaining only terms that will contribute up to single 
logarithmic accuracy and write [20] 

1 /"^sin(6A)e-^°^(^)«e-^°^(i/^) 

TT Jl/A b 



e-ri^R' sec^R' 2^ (-1)" e^^^^+^^T^ 
r(l + R') (2n+l)\2n + l- R' 

(5.5) 



21 



where the function R' is 

i?'(l/A) = ?^ln^. (5.6) 

TT A 

The result above corresponds to a Sudakov behaviour with a peak at i?' = 1. 

However we have only considered the contribution from h > h^. The reason for 
doing this is that attempting to describe the whole h integral by a Sudakov behaviour 
(and retaining only terms up to SL accuracy as above) would produce only the first 
term on the right-hand side of the above equation, which is divergent at i?' = 1. 
This is an indication that near this region, a Sudakov behaviour is not the dominant 
contribution. 

The physical reason for the above statement is mentioned below and appears 
in many other examples including the Drell-Yan pt distribution [27,28]. At very 
low A (below the peak region, R' 1) the dominant mechanism by which one 
obtains a low Et difference between final state jets is vectorial cancellation, rather 
than soft emission. In fact this mechanism is dominant at small A since only a 
one-dimensional cancellation is required, | Xli^^:*! ~^ 0> i'^ present case. Hence 
the Sudakov peak at /?' = 1 is washed out. In problems requiring a two-dimensional 
vectorial cancellation [27], the vectorial cancellation is important only beyond the 
Sudakov peak which therefore appears. 

To see how the recoil cancellation mechanism behaves, we have to consider the 
small h contribution : 



-/ 



^ sin(6A)e-«-('') = - V ^ / dh fc^^e-^B. W (57) 

h TT^ (2n+l)! io ' V ^ 



where we expanded the sine function. In general, this integral has to be done nu- 
merically. However, the simple form of RT)\,{h) allows us to integrate eq. (5.7) term- 
by-term analytically, and obtain 

I *6-e-.'.^i^i-^*(vs4-?^). (5.8) 



where the quantity a and the function ^{x) are 



TT 



/X 
dte-^*\ (5.9) 
-00 



The series above is rapidly convergent. At very small A one obtains the leading 
behaviour 

'^-(A)«(/^e-™eV-j^, (5.10) 

which means that the slope of the resummed dijet distribution near A = goes 
linearly to zero, with a proportionality constant which behaves as which 
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arises from integrating a Sudakov form factor. This behaviour may of course receive 
corrections by terms that start at relative order ag, and we need to account for these 
pieces with a matching procedure. 

The actual value of the coefficient of the linear behaviour at small A will ob- 
viously also depend on the SL terms, but the qualitative behaviour is as we have 
demonstrated at the DL level. 

We next wish to present some plots for the resummed slope cr'{A) for dijet 
production in the Breit frame of DIS. However before doing that we make some 
additional points about the numerical computation. In order to ensure physical 
behaviour over the entire range of b integration, we can redefine the resummation 
variable bQ, by making the replacement — > S((5~"'^a/T+~&^Q^) [29]. This is an 
ad-hoc prescription to some extent but it leaves the soft region (large b behaviour of 
the radiator) unaltered, and does not change the linear behaviour at small A, that 
we described earlier. Other such prescriptions are also possible but would differ only 
in the large A region, which in any case requires matching to fixed order. We also 
point out that to avoid the Landau pole in the running coupling (and unphysical 
behaviour of the parton distributions) we have to put a cut-off on the b integral at 
large b. We take this to be at the Landau pole, A = /^oa^lnfeQ = 0.5, although 
varying the position of this cut near the vicinity of this value makes no noticeable 
difference at all to the results. We shall comment on the role of non-perturbative 
effects in our conclusions. 

In all the presented plots we choose for the DIS variables the values x — 0.01 
and Q — 20 GeV, consider dijet events with the cone algorithm described earher and 
fix the opening angle S — 0.2 with the minimum transverse energy £^min = 10 GeV. 
In order to have two well separated jets which respect the condition S^/[ij] -C 1 , 
where [ij] is the interjet opening angle, we impose a cut |?7| < 1 on the rapidity rj 
of the jets (with respect to the photon direction) in the Breit frame. This works in 
our case, since for the basic Born configuration we integrate over, the jets must have 
equal (cancelling) transverse momenta with respect to the photon direction. On an 
experimental level (or beyond leading order jet production), the cut with respect to 
the photon axis is not a sufficient means to ensure the separation of the high Et 
dijets, and we need to impose a cut on say the interjet rapidity. This is however 
not needed for our purposes here. Additionally we use the MRST2001_1 parton data 
set [30] and the pdf evolution code described in Ref. [23]. 

In figure 3 we plot the ratio 

D{A) = a'{A)/a',{A) , (5.11) 

where crQ(A) is the Born value for the slope (2.6). It is interesting to compare the 
behaviour obtained after resummation (the upper curve), with the logarithmically 
enhanced part of the NLO result, which we have computed. While the resummed 
result has the same sign as the Born term (both are negative as required on physical 
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Figure 3: The ratio D(A) at LL and NLO (soft/collinear) accuracy, for Q = 20 GeV and 
X = 0.01. 

grounds), the NLO result takes the opposite sign to the Born term, at small A, and 
becomes divergent. 

Note also that the resummation predicts a slope that vanishes linearly at small 
A, while the leading (Born) order slope is finite in this limit. The corrections to this 
resummed result will, at small A, at best start at relative order as- Hence we still 
(even after matching to NLO) expect a value of the slope that in the small A region, 
is much smaller than that obtained at leading order. 

Although we do not explicitly indicate it, it should be understood that the in- 
tegration over the hard Born configuration is done with the angular cuts described 
above. The final results include also the gluon incoming channel as well as the 
transverse and longitudinal contributions to the result. 

In figure 4 we can see the impact that NL logarithms have on the slope. There 
we compare the full NLL answer (5.2) with the LL result, obtained from (5.2) by 
neglecting all NLL effects, i.e. setting R{b) = -Rll(&), freezing the scale of the pdf's 
at Q and making the replacement S{b) —>■ 1. One can observe that both curves 
show a linear behaviour for small A but that the LL curve is steeper than the NLL 
one. The non-global component of the NLL result however, makes little difference 
to the final results (actually negligible in the very small A region and less than 5% 
in the whole selected range). The reasons for this were mentioned before and this 
observation is perhaps consoling in light of the fact that these are only known in the 
large limit. 

6. Conclusions 

Here we present the main conclusions from this work. We have demonstrated that 
the region of symmetric Ef cuts in dijet production can be handled with an all order 
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Figure 4: The ratio D{A) at LL and NLL accuracy, for Q = 20 GeV and x = 0.01. 

resummation. The resummation is carried out in the slope of the dijet rate c"'(A). 
The result we obtain shows that the slope of the resummed dijet rate (t'{A) is negative 
throughout and vanishes at small A, at least up to corrections of relative order a^. 
The NLO calculation on the other hand is positive and divergent at small A. 

The resummation we carried out enables us to probe the Et distribution of the 
higher energy jet, arbitrarily close to the cut on its energy, which would otherwise 
be unsafe. Additionally after matching to fixed order, one can in principle achieve 
a better estimate of the rate at A = 0, o'(A) — cr(0) = dxa'{x), where one can 
choose A in the region where the NLO calculation is reliable and hence use the NLO 
value. In fact choosing A large enough that the rate o"(A) vanishes would enable the 
direct determination of cr(0). 

Our calculation was performed using a simple variant of the cone algorithm [15]. 
As we stressed in the paper, other variants of this algorithm can be employed and 
would lead to different results arising from NLL differences. However we still expect 
our result to serve as a qualitative model for the NLL terms, in all cases where the 
recombination scheme ensures that the Et mismatch between the dijets is given by 
the component kx (defined earlier) of soft gluon momenta flowing outside the jets. 
This is because however the jets are defined, soft radiation at large angles to the jets 
will follow a three-particle antenna pattern identical to the one employed here. The 
terms that depend on the geometrical details of the algorithm, i.e. the finite cone 
size pieces, will of course be different but computable at least in the global term. 
For algorithms that involve clustering such as the inclusive kt algorithm with an R 
parameter, the situation with non-global logs will be different. In fact one expects 
the clustering procedure to reduce this component [31], which in any case does not 
significantly affect the result, and hence we should be able to extend our computation 
to apply to that case. 
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We have not mentioned, thus far, that in general one may expect a small Et to 
be generated on the incoming leg, due to non-perturbative effects (intrinsic kt). This 
would also lead to a mismatch between the jet E^s. However this effect is expected 
to correct the radiator at the level of terms quadratic in b, which would lead to a 
power correction, as in the case of the Drell-Yan pt distribution [33]. 

Another important development that is needed is the matching to fixed order 
of our calculation. This would ensure that one can describe the slope at the two- 
loop level completely (accurate even at larger A), if one replaces the pathological 
logarithmic behaviour with our all-order resummation. 

We also mention that it is possible to perform the resummation in other ways, 
to obtain a finite rate with symmetric cuts. For instance one can resum threshold 
logarithms in the invariant mass distribution of the dijet pair. Then one can probe 
the mass distribution safely, even in the highest mass bin, with symmetric cuts. An 
integration over all mass bins would yield a rate that is finite with symmetric cuts. 
This will be the procedure adopted in [32]. 

As further extensions of this work, one can conceive of dijet production via the 
resolved photon contribution rather than DIS. This would involve a four particle 
antenna, rather than one with three particles, as was the case here. Similar issues 
will also arise in the case of dijets produced via 77 collisions at LEP [10] and the 
case of prompt di-photon hadroproduction as we mentioned before [12]. 

Lastly we hope that the work carried out here will eventually lead to comparisons 
with experimental data. However much work remains to be done in terms of matching 
to fixed order, which will include having to consider jet production beyond leading 
order, as well as adjusting the details of the calculation to a different jet algorithm. 
We intend to address these issues in forthcoming work. 
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